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Highlights 
e Layered MgBz2 has the potential to serve as a catalytic and energetic fuel additive 
e This potential can be further enhanced by activating MgB: using ball milling 
e Adding Iwt. % ball-milled MgBo increases the energy of ammonium perchlorate 
by 78% 
e Adding Iwt.% ball-milled MgBz2 reduces the decomposition by ~73°C 
e The first study that utilizes MgB2, a conventional superconductor, as a fuel additive 


Abstract 


In this article, we studied the effect of layered magnesium diboride (LMDB) and its micro and 
nano derivatives as catalytic and energetic additives towards the thermal decomposition of 
ammonium perchlorate (AP). We synthesized micro derivatives by mechanical activation of 
MegBo2 using ball milling and nano derivatives by dissolution recrystallization method. The 


mechanically activated MgB2 (MA-MgBz2) showed better performance than pristine MgB2 (P- 


MgB) and boron-based nanosheets (NSs) on improving the thermal decomposition 
characteristics of AP. Addition of one wt. % of MA-MgB2 to AP remarkably enhances the 
energy release by 78% and significantly reduces the decomposition temperature by ~73°C. 
Using the isoconversional method, we obtained activation energy, pre-exponential factor, and 
reaction model of AP and AP+1%MA-MgBo. The MA-MgB> performs much better than many 
other catalysts at similar loadings. This study provides a newer perspective to utilize the 


untapped potential of metal diborides and their derivatives for composite solid propellants. 
Keywords: layered magnesium diboride, ammonium perchlorate, catalytic additive, and 
composite solid propellants 


1. Introduction 


Ammonium perchlorate (AP) is a widely used energetic material that acts as a strong oxidizer 
and facilitates self-sustained combustions [1]. It is a major ingredient in the currently used 
composite solid propellant fuels (CSPFs) and accounts for ~60-80% of the total weight of 
CSPFs [2-5]. The thermal decomposition characteristics of AP significantly affect the 
combustion performance of the propellant, which in turn determines the degree to which 
energy can be extracted [6]. To improve the thermal decomposition of AP, fuel additives 
and burn rate modifiers are incorporated in the fuel [7]. Generally, metals with a high heat of 
combustion are used as fuel additives [8,9]. Some energetic metal fuel additives include Li, 
Be, B, Si, Al, Mg, Ti, and Zr (heat of combustion values of energetic metals listed in Table 
SI, SD [9,10]. Li and Be are not widely used because they are highly toxic and expensive 
[10,11]. Aluminum is the most commonly used metal fuel additive because of its wide 
availability, ease of handling, the high heat of combustion (31.4 MJ/kg), and also more 


economical [10,12]. 


On the other hand, several burn rate modifiers are added to AP to improve the burning rate 
of the propellant. The burn rate modifiers or combustion modifiers are classified as catalysts 
(increase the burning rate) and inhibitors (decrease the burning rate) based on their activity 
on the combustion performance of the propellant [13]. The burn rate modifiers include 
metals, metal oxides, and mixed metal oxides as they can exhibit variable oxidation states by 
forming complex molecules with the reagents [7,14-16]. The formed metal complex 
molecules can either donate or accept the electrons from the transition state of the reaction 
and facilitate the reaction to occur [14]. Conventionally, transition metal oxides (TMOs) have 
been used as catalysts, and some commonly used catalytic additives include Fe203, CuO, 


MnOdz, Ni2O3, copper chromite, and S102 [7,13,17—19]. LiF is reported as an inhibitor to the 


decomposition reaction of AP, and it decreases the burn rate of the propellant by forming 


NH4gF through the formation of LiF-NH4ClO4 eutectic mixture [13]. 


Boron is another universal fuel additive in propellants; it exhibits the highest theoretical 
heat of combustion in terms of both volumetric (135.8MJ/L), and gravimetric (58.5MJ/kg) 
terms. Its gravimetric heat of combustion is nearly double that of aluminum, and also greater 
than several other traditionally used metal fuel additives [20,21]. Yet, scientists are not able to 
effectively utilize the potential of boron because of the challenges associated with its 
combustion. The challenge is that upon combustion, boron forms a viscous boron oxide layer 
on the surface that limits further oxidation of boron [22,23]. Kuwahara et al. reported that 
increasing the surface area by nanoscaling/micro scaling the boron improves its combustion 
efficiency [24]. However, by reducing the particle size, it is not possible to completely 
overcome the problem of oxide layer passivation. Furthermore, at a small size scale, the 
problem of agglomeration increases. Moreover, the methods to obtain nanoforms of boron are 
not scalable [25-27]. Some studies have reported that the challenges associated with the 
combustion of boron can be overcome by surface modification and coating the individual 
particles with LiF, Viton A, and silane [28-31]. The alternative strategy that was recently 
reported to overcome the challenge is to add active metals to boron. Liu et al. investigated 
various ratios of Mg/B mixtures and shown that the addition of magnesium to boron improves 
the burning of boron particles [32]. Recently, metal borides are gaining renewed attention as 
alternative energetic materials to boron as fuel additives because of their high heat of 
combustion. Several researchers have begun investigating various methods to synthesize high- 


quality metal borides for solid propellant rocket engines [10,33—35]. 


In the past seven years, our research group has established the science of nanoscaling 
metal borides (MBs) of the form MBz2 to obtain boron analogs of graphene. MBs of the form 
MBz2 [36] consists of metal atoms sandwiched between covalently bonded hexagonal 
networks of boron (Figure 1). MBs are commonly used in steel hardening, turbine blades, 
corrosive resistant, and refractory coatings because of their excellent mechanical, chemical, 
and electrical properties [37]. Until our previous studies [38-42], they have never been 
probed as a potential candidate for yielding quasi-2D boron nanostructures. We had revisited 
the chemistry of a well-known layered material, magnesium diboride (MgB2), popularly 
known for its superconducting properties and had shown that it is possible to exfoliate MgB2 
[38—42]. With a curiosity to understand the physico-chemical interaction between MgB2 and 


water, we had also discovered a simple and scalable recipe to obtain various quasi 2D boron- 


based nanostructures in high-yield [40]. The quasi 2D boron-based nanostructures utilized in 
this study were synthesized using the recipe of dissolution-recrystallization method [40] 


discovered by us earlier. 


In this article, we firstly showed the effect of quasi 2D boron-based nanostructures 
on thermal decomposition of AP. We then demonstrated that pristine MgB2 (P-MgBz2) is a 
potential catalytic and energetic additive for enhancing the thermal decomposition 
characteristics of AP. To further enhance the activity of P-MgB2, we prepared mechanically 
activated-MgB2 (MA-MgB> or micro derivatives) by developing an optimized ball milling 
recipe (Fig.1). Finally, we show that the addition of one wt. percentage of MA-MgB>2 
enhances the thermal decomposition of AP leading to an enhanced release of energy, better 
than several other existing catalysts. To the best of our knowledge, there are no existing 
reports on using MgB: and its derivatives as additives for tuning the exothermicity of AP. 
This study provides a maiden insight on utilizing the layered magnesium diboride (LMDB) 
and its micro and nano derivatives as potential catalytic and energetic additives for future 


solid propellants. 
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Figure 1. Schematic representing the formation of mechanically activated-magnesium diboride 
(MA-MgBz2) using planetary ball milling. (a) structure of layered magnesium diboride (LMDB) 
representing the presence of magnesium atoms sandwiched between boron honeycomb planes, 
(b) planetary ball milling of MgB crystals yields in the formation of MA-MgBz, (c) composite 
solid propellant containing AP and one wt. percentage of MA-MgBa», and (d) DSC curve 
showing the effect of one wt. % of boron-based Nanosheets (NSs) and MA-MgB> on the high 
temperature decomposition (HTD) peak of AP. 


2. Materials and Methods 
2.1 Reagents and Chemicals: 


The received ammonium perchlorate (NH4C104) has a mean particle size of ~300 um. Standard 
Magnesium diboride (MgB2) powder of -100-mesh size with => 99% purity was purchased from 
Sigma—Aldrich and used as pristine-MgB2 (P-MgB2). Deionized (DI) water was used for the 
ultrasonication of MgBo. 

2.2 Micronization of Ammonium Perchlorate (AP): 

AP particle size of ~30 um was prepared using a planetary ball mill (Insmart, India), processed 
for 35 min at 375 rpm with a 4:1 ball to powder ratio (zirconia ball diameter is ~10 mm). 
Further, various size particles were segregated by mesh sieving. The particle size analysis of 
micronized AP is shown in Fig. Sl. With this milling parameters, we majorly obtained 30 um 
size and a tiny fraction of 16 um size AP particles. For all our studies, AP30 is only used, and 
until our further studies, the micronized AP particles are sealed tightly and stored under 
vacuum. 

2.3 Synthesis of various boron-based nanostructures: 

Various boron-based nanostructures (NS) were prepared using a sonication process according 
to the dissolution and recrystallization recipe showed in our previous work [40]. In brief, 3 
mg/ml of P-MgB2 in DI water was taken in a double-walled jacketed vessel at 25°C. The 
mixture is probe sonicated using a half-inch probe (Qsonica 500 W) with operating parameters: 
30% amplitude and 10 seconds ON, 10 seconds OFF pulse for 30 minutes. The obtained dark 
suspension was left undisturbed for 24 hours and then vacuum filtered using 0.22um filter 
paper (Durapore GVWP 47 mm). The obtained golden yellow color filtrate was allowed to age 
for Oh, 24h, 48h, and 72h to get boron-based nanodots, nanograins, nanoflakes, and nanosheets 
respectively. The powder form of these boron-based nanostructures was obtained by 
lyophilizing the filtrate at the above specified aging times. 

2.4 Preparation of mechanically activated-MgB2 (MA-MgB2): 


The MA-MgBo were prepared by planetary dry ball milling (Insmart, India), of P-PMgBo. A 
measured quantity of P-MgB> crystals was taken in a zirconia bowl (10 cm inner diameter) to 
which zirconia balls (6mm diameter) were added. The balls to powder ratio were maintained 
as 5:1 (w/w); a milling speed of 150 rpm was used, and the milling time varied from 2h to 30h. 
During milling, at specific time intervals (e.g., 2h, 4h, 6h, 17h, 25h, and 30h), a small amount 


of sample (~1 gm.) was collected and stored in vacuum for further characterizations. The ball- 


milling was performed with different cases to identify the optimal milling parameters. The 
various cases are as follows: 

Case 1: Milling with 6 mm balls 

Case 2: Milling with 10 mm balls 

Case 3: Milling with 90% of 6 mm balls and 10% of 10 mm balls 

Case 4: Milling with 80% of 10 mm balls and 20% of 6 mm balls 

In all the cases, the milling speed is 150 rpm, ball to powder ratio is 5:1 (w/w), and by 
comparing the DSC studies of all milled samples for various cases, we found that milling with 
6mm balls for 17h is effective in terms of energy release and in reducing the decomposition 
temperature. Therefore, all further characterizations were done on the MA-MgBo2 sample 
obtained after 17h of milling from case-1. The particle size analysis of the P-MgB2 and MA- 
MgB: are shown in Fig. S2. With increasing milling time, the average particle of MgB2 was 
found to be significantly decreased. 

2.5 Preparation of composite solid propellant fuels (CSPFs): 

The micronized AP (~30 um) and P-MgB2/MA-MgB2/boron-based NSs were mixed in a 
regular mortar and pestle at a mass ratio of 99:1 to prepare the CSPF samples for thermal 
decomposition analysis. In general, ~1980 mg of micronized AP sample is mixed with ~20 mg 
(one wt. %) of P-PMgB2/MA-MgB2/boron-based NSs in a mortar and pestle and gently mixed 
to obtain a uniform mixture. We also observed that the addition of P-PMgB2/MA-MgB2/boron- 


based NSs improves the free-flowing nature of the micronized AP particles. 
2.6 Activation energy using Kissinger-Akahira-Sunose (KAS) method: 


To determine the activation energy, we followed Kissinger-Akahira-Sunose (KAS) method, a 
model free isoconversional integral method formulated for linear non-isothermal heating 


studies [43,44]. The KAS method is based on an approximation described by equation (1). 





E 
ln Pi = Constant - — (1) 
Poi RT qi 


Where, 

6; =i" heating rate 

T= temperature at which the extent of conversion a is reached under the i™ heating rate 
Eq = activation energy at the extent of conversion a, and 


R = universal gas constant 


A detailed derivation of the equation is shown in the kinetics section of the supporting 
information file. The extent of conversion or reacted fraction a, as a function of time is obtained 
by equation (2). 

m(i) — m(t) 


9° n= MP) 


(2) 
Where m(i) is the initial mass, m(t) 1s the mass at any time t, and m(f) is the final mass obtained 
from TGA analysis. To use this method, we need to obtain TGA data for at least three different 
heating rates. We have performed the TGA analysis at the heating rates (f;) of 5, 10, 15, 20, 
and 25 K/min. By using equation (2), we determined the conversion a(t) for each heating rate 
(Bi) and plotted conversion a(t) vs. temperature (Ty ;) for all the heating rates (f;). We then 
obtained the exact values of temperatures for discrete values of conversion (a = 0.1 to 0.9 with 
a step size of 0.05) by interpolating the data of conversion and temperature at each heating rate. 
For the conversion a = 0.1, we noted the temperatures at the different heating 


rates 61, b2,f3,B,,and Bs as Ty 1,T¢2,Te3,Te¢,4,and Ty 5 respectively. Using the values of [; 


Qa,l 


and Ty ;, we utilized the equation (1) and plotted In (4+) vs. — for the conversion a = 0.1 to 


obtain the value of activation energy (Ea-0.1) from the slope (-E,/R) of the linear fit. By plotting 


1 . , ue 
ln (4) Vs. -— at each discrete value of conversion, we can obtain the values of activation 
al 


energy (Ea) for each conversion. In this way, we obtained the plot of activation energy (Ea) vs. 
conversion (a) for the thermal decomposition of the composite solid propellant AP and 
AP+1%MA-MegB2. This plot is helpful to understand the decomposition behavior at various 


stages of conversion. 


2.7 Master plots based on differential form of generalized kinetic equation 


to determine the reaction model: 


Master plots are theoretical curves that depend on rate determining mechanism and are 
independent of reaction parameters [45]. We utilized the differential form of reaction models 
[f (a) | shown in Table 4, to obtain the theoretical master plot in differential form. The derivation 


for the equation of master plot is shown in supporting information file. 


Generalized equations of differential form of master plot are as shown below: 
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Equation (3) represents that for a given value of a, the experimentally determined value of 


da 


reduced-generalized reaction rate /d9 corroborates with the theoretically calculated 


1/49) os 


AO hen an appropriate differential reaction model [f (a)] is applied. Furthermore, 
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Vs. conversion (a) 


f (0.5) 1s constant for a given reaction model shown in table 4. Plotting a f 


gives the theoretical master plot, which is shown in Figure 8. From equation (4), the reduced- 


da 


reaction rate can be obtained using the experimental data from TGA analysis 


do 
"* 148) anos 
recorded at a constant heating rate. The average activation energy value obtained from any of 
the model free isoconversional methods can be used. Comparing theoretical master plot with 


experimental reduced rate, we found the suitable match of reaction model of the sample. 
2.8 Determination of pre-exponential factor (A) 


After determining the suitable reaction model, we obtained the pre-exponential factor (A) using 
Coats-Redfern method, a well-known integral form of kinetic equation: [43]. The equation is 


as shown below: 











| A,R [. 2RT. E 
In Ca = In(—S_|1 -—?)) (5) 
T2, g(a)E E, |) RT a, 


Where, (6; is the heating rate, and T, ; is the temperature for an i heating rate at a particular 
conversion a. This equation is similar to KAS equation shown in equation (1). So plotting 


1 : . , . 
ln (£:) VS. —— as explained earlier, we obtained the slope and intercept. From the intercept, 


al Qa,l 
using the integral form of the reaction model g(a) determined from the master plot curve, we 


obtained pre-exponential factor (A,,). Here, Tey, 1s the mean temperature calculated from the 


temperature values recorded for different heating rates for a particular value of conversion. In 


this way, we obtained pre-exponential factor at each discrete value of conversion. And the 
constable plot of InA vs. Ea fitting to a straight line shows the kinetic compensation effect at 
every stage of conversion. From the obtained kinetic triplets (E, A, reaction model), we can 


provide the complete reaction rate for the thermal behavior of a solid material. 


Results and discussions 

The quasi 2D boron-based nanostructures were obtained by dissolution and 
recrystallization of MgBo crystals (Fig. la) in water using the protocol described in our 
previous work (brief about synthesis in the methods section) [40]. The morphology of the P- 
MgB: and the quasi 2D boron-based nanostructures (nanodots (NDs), nanograins (NGs), and 
nanosheets (NSs)) are shown in Fig. S3. To see the effect of these boron-based nanostructures 
on thermal decomposition of AP, we have performed the simultaneous thermal analysis 
(preparation of AP and its composites (CSPFs) in the methods section). Fig. 2a shows the 
differential scanning calorimetry (DSC) graph of pure AP or with no additive, which typically 
indicates one endothermic peak at 245 °C corresponding to the solid-state phase transition from 
orthorhombic to cubic phase [2]. The other two exothermic peaks in the range of 270 to 440 
°C correspond to the decomposition peaks of AP. The small exothermic peak at 283 °C 
corresponds to low-temperature decomposition (LTD) resulting from the partial decomposition 
of AP, and the higher exothermic peak at 428 °C, which is attributed to the complete 
decomposition, is known as high-temperature decomposition (HTD) peak [46]. In general, the 
thermal decomposition reaction of AP is described as follows [2]: 


NH,ClO, @ NHz + ClO, @& NH3g) + HClO4(g) 
= N H3/s) + HClO, s) (6) 


The two-stage thermal decomposition of pure AP was also verified through 
thermogravimetric analysis (TGA). Fig. 2b shows the mass loss of pure AP; at LTD and HTD, 
the mass loss corresponds to 24% and 75%, respectively. The detailed information of AP 
decomposition data is summarized in Table 1. Fig. 2a shows the DSC curves of AP 
decomposition in the presence of one wt. percentage of boron-based nanosheets (NSs); Fig. S4 
shows these curves in the presence of NDs and NGs. With the addition of one wt. percentage 
of boron-based nanostructures (1.e., either NDs/NGs/NSs) to AP, we observed a right shift of 
~20 °C in the LTD peak of AP and no change in either phase transition peak or HTD peak as 
per DSC analysis. This signifies that the initiation of AP decomposition is delayed with the 


addition of boron-based nanostructures. This is probably due to the presence of functional 


groups on the nanostructures [40], which is evident from the FTIR data shown in Fig. S5. These 
boron-based nanostructures were found to be functionalized with chemical functions such as 
hydrides, oxides, and hydroxide functional groups, which might be hindering the LTD process 
of AP [2]. This is also evident from the TGA mass loss curves as we can observe a slight 
deviation at both the stages of mass loss as compared to pure AP (Fig. 2b and S4). We observed 
clear shifts in the decomposition peaks of AP and its composites in the DTG plots. The effect 
of NSs as additives on thermal decomposition of AP is shown in Fig 2, and the effect of NDs 
and NGs are shown in Figure S4. We observed that LTD peak is right shifted in the presence 
of nanostructures. However, while the position of HTD peak for AP does not change in 
presence of NDs and NGs but it is significantly left shifted in presence of NSs. This suggests 
that the presence of boron-based nanostructures hinders only LTD and promotes HTD, which 
is clearly evident from the DTG plot shown in Fig. 2d in the presence of NSs. We also observe 
an increase in the exothermicity of AP from 763 J/g to 1094 J/g (Fig. 3a and Table 1), 1.e., a 
~4()% increase in the energy release of AP in the presence of boron-based nanosheets. To verify 
this, we have performed the hyphenated TGA-mass spectrometry (MS) of gas residues 
produced that will demonstrate later the boron gaseous species released during the HTD 
process. 

To obtain more insights into the performance of these boron-based nanostructures derived from 
LMDB, we also conducted control experiments with the parent material, P-MgB2 as an additive 
in AP decomposition. Surprisingly, we found a significant enhancement in the decomposition 
characteristics of AP in the presence of P-MgBo. We noticed a sharp rise in the HTD peak of 
AP along with a shift towards lower temperature. Almost no changes in the phase transition 
and LTD peaks were observed in the case of P-MgB2. We observed that the HTD of AP is 
significantly promoted, just by an addition of one wt. % of P-MgBo (Fig. 2a). We noticed a left 
shift of ~62°C in HTD peak of AP along with an increase in the exothermicity of AP from 763 
J/g to 1239 J/g (Fig. 3a and Table 1), 1.e., a ~62% increase in the energy release of AP in the 
presence of P-MgB>o. This unprecedented increase in exothermicity and decrease in 
decomposition temperature of AP can be attributed to the combined effect of presence of 
magnesium metal atoms (as catalyst) and boron planes (as fuel additive) together in the P- 
MegBz2. This indicates that P-MgBz2 acts as a catalytic as well as energetic additive towards the 
thermal decomposition of AP. To the best of our knowledge, there are no existing reports, 
which study the effect of MgBz2 as an additive on thermal decomposition characteristics of AP. 
Motivated by the superior performance exhibited by MgBz, we endeavored to develop a way 


to nanoscale MgB2 with minimal or no functionalization, as nanomaterials are expected to 


show a superior activity than their parent material [47]. Based on the knowledge obtained from 
our previous studies on exfoliating MgB2 [38-42], we know that liquid phase exfoliation yields 
in the functionalized form of boron-based nanosheets. So, here we decided to use the dry ball 
milling method [48] to nanoscale MgB: as this method 1s easily scalable and is used to obtain 
pristine or minimally functionalized nanosheets from other layered materials. 
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Figure 2. (a) DSC curves, (b) TGA curves, and (c), (d) DTG curves for the LTD and HTD of 
AP decomposition in the absence and presence of additives derived from LMDB at a 1% mass 
basis — No additive or Pure AP, Nanosheets (NSs), Pristine (P)-MgB2, Mechanically 
Activated (MA)-MgB2-2h, MA-MgB2-4h, MA-MgBo-6h, MA-MgBo-17h, MA-MgB2-25h, 
and MA-MgB2-30h, respectively. Among all, the presence of MA-MgBo-17h significantly 
increased the intensity of the HTD peak of AP along with increased energy release and 
decreased decomposition temperature. 

In this pursuit, we identified an optimized set of milling parameters such as milling 
speed, milling time, ball diameter, and ball to powder ratio. During milling, the balls have an 
impact of compressive and shear forces on the milling material. If compressive forces 
dominate, then it results in more fragmentation of the material while shear force dominant 


milling results in exfoliation of layered material. Li et al. have shown that milling of bulk h- 


BN at low speeds results in shear force dominant milling and produces a high yield of h-BN 
nanosheets [49]. Considering these facts, we designed four cases of ball milling to identify the 
optimum milling parameters, such as ball diameter and milling time. We fixed the milling 
speed as 150 rpm (low milling speed for shear force dominance) and ball to powder ratio as 
5:1 (w/w) and performed the milling of P-MgBz for all the four cases. The milled samples were 
collected at different time intervals (O, 2, 4, 6, 17, 25, and 30h) for all the four cases and 
subsequently used to prepare the composites of AP. We performed the DSC analysis of all the 
composites (Tables S2-S5). We found that milling with 6 mm balls has shown better 
performance and we also have observed a gradual improvement in performance with increasing 
the milling time, as shown in Fig 2. Therefore, further analysis was done on the milled MgBo 
samples obtained from milling with 6 mm balls. From here, we termed the milled samples as 


mechanically activated-MgB2 (MA-MgBz). 
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Figure 3. (a) Energy released and (b) HTD peak temperature of the AP decomposition in the 
absence and presence of additives derived from LMDB at a 1% mass basis — No additive or 
Pure AP, Nanosheets (NSs), Pristine (P)-MgB2, Mechanically Activated (MA)-MgB2-2h, MA- 
MgBo-4h, MA-MgBo-6h, MA-MgBo-17h, MA-MgBo-25h, and MA-MgBo-30h, respectively. 
In the presence of MA-MgB>-17h, the energy released 1s the highest with 1355 J/g along with 


a significant reduction in HTD peak temperature. 


Table 1. Results of thermal decomposition of AP in the absence and presence of one wt. % of 


additives derived from LMDB at a heating rate of 10 K/min 


Peak End The energy 
Temperature decomposition Mass loss released 


Composite Tm/°C  temperature/°C Am/% AH/J g" 


LTD HTD Am1 Am2 


No additive 283 428 435 24 i 763 
F-B-NDs 305 431 433 28 71 1050 
F-B-NGs 306 429 432 26 73 1061 
F-B-NSs 304 430 433 ae dz 1094 
P-MgBo 278 366 37] - 100 1239 
MA-MgB2-2h 276 366 371 - 96 1158 
MA-MgB2-4h 276 365 370 - 96 1251 
MA-MgB>-6h 292 370 aS - 95 122] 
MA-MgB2-17h ZA] 362 366 - 93 1355 
MA-MgB2-25h 288 355 360 - 98 L221 
MA-MgB2-30h 286 361 367 - 96 1274 


To investigate the degree of exfoliation upon milling, we collected the powder X-ray 
diffraction (XRD) patterns of the samples obtained at different milling times (Fig. 4a). The 
XRD spectra of the P-MgB2 crystals and the 2 h, 4h, 6h, 17 h, 25 h, and 30 h milled samples 
are compared in Fig. 4a to see the changes in peak intensities and positions. It can be seen that 
with an increase in the milling time, the peaks exhibit a slight reduction in the intensity and no 
shift in the peak positions were observed. In a crystal structure, the multiples of (O01) plane 
correspond to the stacking of layers in c-direction, and the (hkO) planes represent the order 
within the basal sheets [50]. Therefore, to view the effect of milling on c-direction planes of 
LMDB, we compared the enlarged XRD spectra of (O01) and (002) peaks corresponding to the 
P-MgB2 and MA-MgBo-17h (Fig. 4b). It can be seen that even after 17 h of milling, no 
significant reduction in the intensity of both (O01) and (002) peaks was observed. This is 
indicating that exfoliation has not occurred and resulted in the formation of MA-MgB> crystals. 
Extending the milling time further has an inverse effect of milling, which can be evident from 
the slight increase in the intensity of (001) peak for 25h and 30h milled samples in comparison 
with 17h milled sample (Fig. S6). 


b) iim P-MgB, 
—=— MAMgB,-17h 


101) 


(002) (110) 





(002) 


MA -MogB.-17h 


(001) 


Intensity (a.u.) 
Intensity (a.u) 


MA-Mgf.,-2h 








20 30 70 80 90 25 2 £=51  -§3 
20 (degree) 


50 60 
20 (degree) 

Figure 4. XRD analysis of (a) Pristine (P)-MgB2 compared with Mechanically Activated 

(MA)-MgB2 obtained at different milling times; (b) the enlarged spectra of (O01) and (002) 


peaks corresponding to P-MgB2 and the MA-MgB>-17h. With increasing the milling time, the 
P-MgBo2 did not result in significant exfoliation but yielded in the mechanical activation of 
MgB2. 

Therefore, we acknowledge the fact that dry ball milling did not yield in exfoliation. 
However, the MA-MgB> as an additive shows superior activity than P-MgB2 in enhancing the 
thermal decomposition of AP. Fig. 2a shows the decomposition of AP in the presence of one 
wt. percentage of MA-MgB2 obtained at various milling times. We found that MA-MgB2 
obtained by 17 h of milling, enhances the exothermicity of AP in an excellent way better than 
many other existing catalysts (comparison with other catalytic additives is shown in Table 2). 
This remarkable catalytic activity is due to the presence of high-energy material boron (present 
along with the metal atom Mg), combined with the fact that this construct is available in a 
layered form that allows enhanced access to the active sites. This is of immense prospect for 
utilizing in the practical applications of missile fuels. We also determined the net amount of 
heat released in each case, and it was found to be 763, 1239, and 1355 J/g for pure AP and 
composites of AP with P-MgB2 and MA-MgB>-17h samples respectively. The energy released 
data for the remaining samples are tabulated in table |. This shows that we can enhance the 
decomposition energy of AP by ~78% just by adding one wt. percentage of MA-MgB2-17h. In 
addition to the increase of net heat released, we are also able to lower the temperature of 
complete decomposition peak (HTD peak) from 428 °C to 362 °C, a decrease of 66°C (Fig. 3b 
and Table 1). This suggests that MA-MgB>-17h acts as a catalytic and energetic additive in the 
thermal decomposition of AP. Hence, further studies were performed on MA-MgB2 obtained 


at 17 h of milling, and from here, we term it as MA-MgB2. 


Table 2. Comparison of catalytic and energetic activity for AP composites in the presence of 


various additives that are reported in the literature 


Enhancement The shift in 

Loadin in energy HTD peak 
Additive g/wt. % release/% —temperature/°C Reference 
P-MgB2 1 +62 -62 This work 
MA-MgB2-17h 1 +78 -66 This work 
MA-MegB2-25h 1 +60 -73 This work 
Cu@h-BN 2 +36 -45 Huang et al.[51] 
GO D) - -32 Memon et al.[52] 
Al (25 nm) 5 -3] -8 Liu et al.[53] 
Anatase T102 ] +69 -30 Reid et al.[54] 
T102 (20 nm) ] -10 +6 Dixit et al.[55] 
Fe203 (< 40 nm) ] +18 +6 Dixit et al.[55] 
S102 | +5 -2 Dixit et al.[55] 
CuO-h-BN 3 - -50 Paulose et al.[56] 


The decomposition of AP in the presence of P-MgB2 and MA-MgB: as a catalytic 
additive was further investigated by thermogravimetric analysis (TGA), as shown in Fig. 2b 
and Table 1. With the addition of one wt. percentage of P-MgB2 or MA-MgB2 milled samples 
to AP, we observe a significant reduction in the end decomposition temperature of AP. The 
shift in the LTD and HTD peaks can be noticed clearly from the DTG plots shown in Fig. 2c 
and 2d. The MA-MgBo sample shows remarkable catalytic activity by shifting the end 
decomposition temperature of AP from 435 °C to 366 °C (see Table |). Compared with pure 
AP, there is a steeper decline in the TGA curves of samples with P-MgB2 as well as with MA- 
MgB2 samples. This again suggests that MgB> and its derivatives are excellent catalytic 


additives for thermal decomposition of AP. 


To re-confirm the catalytic activity of MA-MgB> on the thermal decomposition reaction 
of AP, we need to know the kinetic triplets (activation energy, pre-exponential factor, and 
reaction model) for the thermal decomposition reaction of AP in the presence and absence of 
additives. So we have performed TGA-DSC studies of the AP and AP/MA-MgB2 composites 
at different heating rates of 5, 10, 15, 20, and 25 K/min. It is known that the thermal 


decomposition of solid propellants is a multi-step process, and the reaction mechanism changes 


with variation in temperature, and thus the activation energy of the system varies [57]. The 
temperature of decomposition and the heat release are highly dependent on the heating rate 
[58]. Fig. 5a and 5b, shows the DSC curves of pure AP and AP+1%MA-MgB> at different 
heating rates, indicating the shift in decomposition peaks to higher temperatures with increase 
in heating rate. All these measurements were recorded in an aluminum pan closed with a 
pierced lid under nitrogen as purge gas with a flow rate of 80 ml/min. From the TGA data 
obtained at each heating rate, we have calculated the extent of conversion (a) using equation 
(2). Fig. 5c and 5d shows the conversion (a) vs. temperature plot of pure AP and AP+1% MA- 
MgB: at heating rates of 5, 10, 15, 20, and 25 K/min. The curves of different heating rates are 


parallel to each other and indicating the delay in conversion with increase in heating rate. 
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Figure 5. DSC curves and Conversion (a) vs. temperature plots for the thermal decomposition 
of (a), (c) pure AP; and (b) (d) AP+1% Mechanically Activated (MA)-MgBo. These plots were 
obtained at different heating rates: 5, 10, 15, 20, and 25 K/min in aluminum pan closed with a 


pierced lid; the measurements were made in a nitrogen environment (flow rate of 80 ml/min). 


To understand the effect of MA-MgB2, we have shown the plots of a vs. T for AP 
decomposition with and without additives at each heating rate in Fig. 6. From the plots shown 
in Fig. 6, we can clearly distinguish that in the presence of MA-MgBz, the curve is highly steep 
and the complete conversion is attained quickly at lower temperatures compared to pure AP. 
We provide a comparison (table 3) to see the difference in conversion at different temperature 
ranges for pure AP and AP+1% MA-MgBp at 5 K/min. From Table 3, we can find that at 
temperatures < 330°C, the conversion is very low and pure AP performs better than the 
presence of MA-MgB>. But from 330-380°C, the conversion spikes up suddenly from 0.19 to 
0.95 in the presence of MA-MgB> whereas pure AP exhibits a conversion only from 0.3 to 0.5. 
This indicates that while pure AP decomposes 50%, the AP with one wt. % MA-MgBz achieves 
95% decomposition, which indicates the merit of using MA-MgB> as an additive to CSPFs like 
AP. Decomposition of last 5% of the material is very slow and requires high temperatures; this 


additive appears to have no effect in in the conversion range of 0.96-0.99. 
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Figure 6. Conversion (a) vs. temperature plots for pure AP and AP+1%MA-MgB: at different 
heating rates: (a) 5 K/min, (b) 10 K/min, (c) 15 K/min, (d) 20 k/min, and (e) 25 K/min. The 
measurements were made in an aluminum pan closed with pierced lid under a nitrogen 


environment (flow rate of 80 ml/min). 


Table 3. Comparison of conversion of thermal decomposition of AP in the absence and 


presence of 1% MA-MgB: at 5 K/min 


Temperature range (°C) 


Conversion (a) 


PureAP AP+1% MA-MgB2 

30-250 < 0.05 < 0.02 
250-280  0.05-0.14 0.02-0.15 
280-330 —0.14-0.30 0.15-0.19 
330-380 —-0.30-0.50 0.19-0.95 
380-412 —-0.50-0.96 0.95-0.96 
412-480 — 0.96-0.99 0.96-0.99 


Now by using the conversion (a) vs. temperature data of different heating rates, we can find 
the first kinetic triplet—activation energy (£,) at various stages of conversion. E,, at each 
conversion is obtained using an integral model free isoconversional method, Kissinger- 
Akahira-Sunose (KAS) method (the data processing is explained in the section 2.6). Following 
the KAS method, we obtained the activation energies of decomposition process with extent of 
conversion for the pure AP and AP+1%MA-MgBz and the respective curves with error bars 
are shown in Figure 7. 
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Figure 7. Activation energy (E£,,) at discrete values of conversion obtained by Kissinger- 
Akahira-Sunose (KAS) method, an integral model free isoconversional approach. The catalytic 
behavior of MA-MgB2 on thermal decomposition of AP is observed at lower (a < 0.2) and 


higher (a > 0.6) conversions. The activation energy curves of pure AP and AP+1% MA-MgB2 


follow different paths with increase in conversion, suggesting completely different thermal 


behaviors. 


From Fig. 7, it is evident that the decomposition of AP and its composite is not a single step 
process as the activation energy is varying throughout the process of decomposition. 
Comparison of Eq vs. a curves for pure AP and AP+1% MA-MgBz, suggests that the thermal 
decomposition in both the cases follows entirely different paths indicating distinct thermal 
behaviors. For the case of pure AP, the activation energy varies from ~78 to 248 kJ/mol 
whereas for the case of AP+1% MA-MgBz», the activation energy varies from ~53 to 212 
kJ/mol. Activation energy values are highly dependent on particle size, heating rate, purge gas 
flow rate, and open/closed lid while performing the decomposition studies of AP. So direct 
comparison with the available literature is not appropriate. Mallick et al. have recently studied 
the thermal decomposition kinetics of pure AP and also provided a comparative plot of 
activation energy vs. conversion of AP from various studies [6]. Among all the studies, the 
decomposition behavior of AP (Ea vs. a) observed by Zhu et al. [59] looks similar to the 
decomposition behavior of our AP+1%MA-MgBz. And the decomposition behavior of AP (Eq 
vs. a) observed by Vyazovkin looks similar to our case of pure AP [60,61]. Overall, we can 
comment that at lower (a < 0.2) and higher (a > 0.6) conversions, the catalytic behavior of 
MA-MgB2 on thermal decomposition of AP is observed. Further, to determine the second 
kinetic triplet-pre-exponential factor (A), we need to first find the reaction model. The reaction 
model describes the rate determining mechanism of the reaction. Using the master plots we can 


find the reaction model (details explained in section 2.7). Fig. 8 shows the theoretical master 
f(a) 





vs. conversion (a) of various solid-state reaction mechanisms listed in Table 4. 


Table 4. Reaction models of most common reaction mechanisms in solid-state reactions 


[6,45,61—63] 
S. : Model 
No Reaction Model aaa f(a) g(a) 
1 Power law p2 a2 a 
2 Power law P3 203 a 
3 Power law P4 haz an 
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Figure 8. Theoretical master plot in differential form representing f(a)/f (0.5) vs. conversion 


(a) for different solid-state reaction kinetic models presented in Table 4. 


To find the suitable rate determining mechanism for our samples, in Fig. 9 we compare the 


da 
theoretical master plots of f(@)/f(0.5)with experimental master plots of /d9 


7 
d 
“/46) «0s 


function of o according to equation (4). From Fig. 9a, the comparative plot for pure AP suggests 
that the experimental master plot is in good agreement with the D3 reaction model implying 
that the rate determining mechanism is three-dimensional (3D) diffusion. However, for the 
AP+1% MA-MgB2 sample, the data for which is shown in Fig. 9b, the comparative plots 
suggest that the experimental master plot follows A2 reaction model implying that the rate 
determining mechanism is two-dimensional (2D) nucleation. This suggests that with the 
addition of one wt. % of MA-MgB> to AP has altered the reaction mechanism by shifting the 


rate determining step from 3D diffusion to 2D nucleation. 
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Figure 9. Comparison of theoretical master plots with experimental reduced rate data using a 
generalized kinetic equation of differential form (Eq.4) for: (a) Pure AP, and (b) AP+1% MA- 
MegB2. The experimental reduced rate data are calculated at a constant heating rate of 10 K/min 
for pure AP and 5 K/min for AP+1% MA-MgBz by using the average activation energy values 
obtained from KAS method. We found that pure AP follows D3 reaction model (Three- 
Dimensional Diffusion) and AP+1%MA-MgB3z follows A2 reaction model (Two-Dimensional 


Nucleation). 


By using the observed rate determining reaction mechanism model from master plots and the 
activation energy values obtained from the KAS method, we have determined the pre- 
exponential factor (A) at discrete values of conversion using the well-known Coats-Redfern 
method [43]. The details of the data processing are explained in section 2.8. The values of (A) 
for pure AP varies from ~7.50x10° to 9.56x10"’ min”! and for AP+1% MA-MgB> it varies from 
~4.67x10"° to 9.95x10"° min’. The linear relationship between the logarithmic form of pre- 


exponential factor (A) and the activation energy (E) for pure AP and AP+1% MA-MgB2 


samples are shown in Fig. S9. The linear regression fit shows the value of R*=0.999 for both 
the samples indicating that for all the heating rates the value of Ag is in accordance with the 
activation energy values at all stages of decomposition process. The range and the average 
values of activation energy(E,), pre-exponential factor (A), and the rate determining 


mechanism for AP and AP+1% MA-MgBz samples are shown in Table 5. 


Table 5. Kinetic triplets (kinetic parameters (Ea, Aa) and reaction model) for the thermal 
decomposition of pure AP and AP with one wt. percentage of MA-MgB2 composite. Activation 
energy values obtained using the KAS method and pre-exponential factor obtained using 
Coats-Redfern method by fitting the reaction models (shown in Table 4) using differential 


master plots. 


Reaction 

Additive E./kJ mol! A/min'! 
model 

No Range: 7.50x 10°+0.52x 107 — 
Range: 78.62 +6.63 — 248.30+24.77 D3 

additive 9.56x10!"46.45x10° 

MA- Range: 4.67x10!943.81x10° — 
Range: 53.0746.26 — 212.82+8.84 A2 

MgB 9.95x10!°+2.06x 10° 


To obtain further insights, we also investigated the decomposition products of AP in 
the presence of boron-based nanosheets, P-MgB2, and MA-MgB> using thermogravimetric 
analysis coupled with mass spectroscopy. Fig. 10a shows the mass loss curve along with the 
gaseous species evolved during the decomposition of AP. The pure AP on thermal 
decomposition was found to release H, NH, OH, H20, NH3, NH4, NO, Cl, HCl, N2O, NOs, 
ClO, HOCI, ClO2, HCIO2, Cle, ClO3, HC1O3, and HClO, species during LTD and HTD, these 
matches well with the reported literature [2,6,64]. Mallick et al. did not observe the detection 
of chloric and perchloric acids and the higher chlorine oxides where as Heath et al. have shown 
the detection of almost all the species similar to our case shown in Table 6. Only the prominent 
species are shown in the mass spectra plot in Fig. 10 and all the gaseous species that were 
detected are listed in Table 6. The NH4 or H20 species curve is not shown in the graph as it is 
overlapping with the curve of OH or NH3. With the addition of one wt. percentage of boron- 
based nanosheets to AP (Fig. 10b), we observed the release of additional gases B3H2, BO2 or 
B3Hi0, and (BO)2 or B4Hio during the decomposition phases of AP. This suggests that the 


boron-based nanosheets releasing B3H2, BO2 or B3Hi0, and (BO)2 or B4aHio species facilitate 
the HTD reactions of AP, which was also evident from the TG-DSC analysis shown in Fig. 2, 
and the additional release of gases is contributing to the increased energy release (as shown in 
Fig. 3a). 
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Figure 10. TGA-MS signals of evolved gases upon thermal decomposition of AP in the 
absence and presence of one wt.% of additives derived from LMDB at a heating rate of 10°C 
min! — (a) Pure AP or No additive, (b) Nanosheets (NSs), (c) Pristine (P)-MgB2, and (d) 
Mechanically Activated (MA)-MgB2 respectively. In the presence of MA-MgBz, a large 
number of species were strongly detected indicating, that it favors the decomposition and 
presence of more active sites enhances the decomposition of AP by promoting the charge 


transfer or redox reactions of AP at HTD by acting as an electrically conductive surface. 


In the case of P-MgB2/MA-MgB> as an additive to AP, we also observed the release of 
additional gaseous species of boranes as well as boron oxides at both LTD and HTD (Fig. 10c 


and 10d). This indicates that P-MgB2/MA-MgB: are not only helping in promoting the charge 


transfer or redox reactions of AP during HTD by acting as an electrically conductive surface, 
but they are also participating in the reaction, which is resulting in the release of boron oxide 
and boron hydride gases. The evolution of BO2 gas was also observed by Liu et al. while testing 
the effect of magnesium on the burning of boron particles [32]. This indicates that P- 
MgB2/MA-MgB2 acts as both catalytic and energetic additive in promoting the thermal 
decomposition of AP, which is also evident from the TG-DSC analysis shown in Fig. 2. The 
relative intensity of gaseous species detected in the case of MA-MgB> is comparatively higher 
than the case of P-MgBo. This promoting effect is attributed to the presence of more active sites 
or an increase in surface area of the MA-MgBz», which was also confirmed through the BET 
specific surface area analysis (shown in Fig. S8 and Table S6). We found that these are 
mesoporous in nature and the specific surface area of MA-MgB> is found to be double that of 
P-MgBz2 and also has increased pore volume compared to P-MgBz. The effect of the same can 
also be witnessed from the sharp rise in the HTD peak of AP in the presence of MA-MgB2. 
The profound increment in energy release along with the completion of decomposition at 
significantly lower temperatures indicates that MA-MgB> 1s a superior catalytic and energetic 
additive. 

Table 6. Comparison of evolved species detected in TGA-MS during thermal decomposition 


of AP in the absence and presence of one wt. % of additives derived from LMDB 


m/z Species No NSs P-MgB2 MA-MgB2 
Additive 

1H v v v v 
15 NH v v v v 
17 OHor NH3 v v v v 
18 H20 or NH4 v v v v 
30 NO v v v v 
35 Clor B3Ho v v v v 
36 HCl v v v v 
43 BO> or B3Hi0 — v v v 
44 NO v v v v 
46 NO» v v v v 
51 ClO Vv v v v 
52 HOCI v v v v 
54 (BO): or B4aHio — v v v 


62 H3BO3 = v v v 
67 ClO2 v v v v 
68 HClO2 v v v v 
70 Ch v v v v 
83 ClO3 v v v v 
84 HClO; v v v v 
100 HClO, v v v v 


The mechanism of AP decomposition in the presence and absence of MA-MgB> as a catalytic 
and energetic additive is explained through a schematic shown in Fig. $10. The presence of 
MA-MgB>z does not affect the phase transition temperature and LTD, its effect is pronounced 
during the HTD reaction of AP. The decomposition mechanism of AP is not clear even today; 
however, the mechanism proposed by Boldyrev is widely accepted. According to Boldyrev, 
the initiation of decomposition of AP is through proton transfer mechanism from cation NH4* 
to anion ClOg, which form NH3 and HCIO4. These intermediate products are partly adsorbed 
on the surface of AP, which undergo further decomposition and oxidation reactions to form 
various gaseous species. While some of the intermediate products directly sublimate without 
undergoing further decomposition and oxidation reactions, the intermediate products pile up 
on the surface of AP at low temperatures and hinder the decomposition reaction, which is why 
we observe a partial low temperature decomposition peak of AP. At higher temperatures, this 
closed gaseous layer bursts, which then leads to the complete decomposition of AP. In the case 
of pure AP, the LTD occurs below 350 °C and HTD occurs above 410 °C and forms various 
gaseous species, as shown in Fig. S10. In the presence of MA-MgB> on the surface of AP, we 
believe that the MA-MgBz2 acts as a catalytic substrate and helps in the redox reactions by 
facilitating a charge transfer. This could lead to decomposition at temperatures below 365 °C. 
This significant shift in the decomposition temperature of AP in the presence of MA-MgB> 1s 
likely due to the extremely conducting nature of MgB2. We also acknowledge that the MA- 
MgB: participates in the reaction forming boron oxides and boron hydride gaseous species, 
which adds more energy to the overall decomposition reaction of AP. Both these aspects make 


the MA-MgB>2 a superlative catalytic and energetic additive for fuels. 
3. Conclusion 
We have studied the effect of layered MgB2 and its derivatives on the thermal decomposition 


characteristics of AP. The nano derivatives were synthesized by a dissolution recrystallization 


method while micro derivatives were obtained by low energy dry ball milling method. With 
the addition of one wt. percentage of these derivatives to AP, we observed promoting effect in 
the decomposition of AP. Among boron-based NSs, P-MgB2, and MA-MgBz, the performance 
of catalytic and energetic activity is as follows: MA-MgB2>P-MgB2>boron-based NSs. With 
the enhanced surface area and mesoporous nature of MA-MgBz, we observed an increase in 
the energy release of AP by 78%, a decrease in the decomposition temperature by ~73°C.The 
variation of activation energy throughout the decomposition reaction is obtained by following 
integral model free isoconversional method. The results suggest that the MA-MgB> is a 
promising candidate producing more effect than the nanostructure catalysts fabricated by a 
sonication process and can be used as a future potential catalytic and energetic additive in 
tuning the exothermicity of AP. Synthesizing pristine or minimally functionalized MgBo 
nanosheets and utilizing them as additives would further enhance the exothermicity of AP and 
it would be interesting to see their effect on decomposition characteristics of AP. This study 
provides a fundamentally new perspective that MgB2 and its derivatives can be used as catalytic 
and energetic additives in the field of solid propellants. It would be promising to explore the 
effect of other members of the MgB2 family (for example, AlB2 and T1B2) on the exothermicity 
of solid propellants. 
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